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I. INTRODUCTION 



Weinreb and Crosby (1973), (henceforth abbreviated WC), 
presented a useful technique for estimating atmospheric 
moisture profiles to assist in obtaining more accurate tem- 
perature retrievals from satellite soundings. Their estima- 
tion procedure was designed to derive the mixing ratios W at 
the lowest 40 levels of the 100-level atmosphere (McMillin, 
et al, 1973) defined in terms of linear j-increments of P 2 ^ 
given by 



The saturated mixing ratio values W g (k) at 11 key levels 
are chosen as predictors from among the lowest 40 j-ievels 



.01,..., 1000 mb in the manner shown in Table 1. Table 1 
also indicates the 11 fixed k-levels designated for computa- 
tional use of the W (k)-values employed as predictors. The 

s * 

data set used by WC to establish the regression coefficients 
was deduced from 1100 island and coastal radiosondes with 
humidities reported to 200 mb. The WC radiosonde sample was 
distributed in an approximately random manner with respect 
to season and latitude. 

The forty regression equations of WC were presented in 
the form 



P . 

3 



. 01(1. +( . 26087836) (j-1) 7 / 2 



( 1 ) 



of eq. (1). In (1), as j = 01 



100, P. ranges over 
J 




y • • • > 



40. 



( 2 ) 
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Table 1. The 100 level atmosphere. Each j is an increment of P 2 / 7 , 
j = 1,...100. Asterisked values are used for predictor 
levels in the regression analyses conducted. 



Level 

j 


Pressure 

p(j) 


Level 

j 


Pressure 

p(j) 


Level 

j 


Pressure 

pa) 




(mb) 




(mb) 




(mb) 


1 


.0100 


35 


30.2057 


68 


271.2453 


2 


.0225 


36 


33.0936 


69 


284.8862 


3 


.0434 


37 


36.1735 


70* 


299.0103 


4 


.0756 


38 


39.4530 


71 


313.6275 


5 


.1219 


39 


42.9395 


72 


328.7482 


6 


.1857 


40 


46.6407 


73 


344.3825 


7 


.2703 


41 


50.5643 


74 


360.5408 


8 


.3794 


42 


54.7182 


75 


377.2336 


9 


.5168 


43 


59.1104 


76 


394.4712 


10 


.6866 


44 


63.7487 


77* 


412.2642 


11 


.8928 


45 


68.6414 


78 


430.6232 


12 


1.1398 


46 


73.7966 


79 


449.5589 


13 


1.4322 


47 


79.2225 


80 


469.0820 


14 


1.7746 


48 


84.9277 


81 


489.2033 


15 


2.1718 


49 


90.9204 


82* 


509.9337 


16 


2.6288 


50 


97.2092 


83 


531.2841 


17 


3.1507 


51 


103.8027 


84 


553.2655 


18 


3.7426 


52 


110.7097 


85 


575.8889 


19 


4.4100 


53 


117.9389 


86 


599.1655 


20 


5.1584 


54 


125.4991 


87 


623.1065 


21 


5.9933 


55 


133.3993 


88 


647.7231 


22 


6. 9206 


56 


141.6484 


89 


673.0266 


23 


7.9461 


57 


150.2557 


90* 


699.0285 


24 


9.0758 


58 


159.2302 


91 


725.7401 


25 


10.3158 


59 


168.5812 


92* 


753.1729 


26 


11.6723 


60 


178.3181 


93 


781.3385 


27 


13.1517 


61 


188.4501 


94* 


810.2485 


28 


14.7804 


62 


198.9868 


95 


839.9147 


29 


16.5049 


63* 


209.9378 


96* 


870.3486 


30 


18.3920 


64 


221.3126 


97 


901.5623 


31 


20.4283 


65 


233.1209 


98* 


933.5674 


32 


22.6209 


66* 


245.3726 


99 


966.3760 


33 


24.9765 


67 


258.0774 


100* 


1000.0000 


34 


27.5024 
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Equation (2) gives the format of the prediction equation 

where j is the level of the mixing ratio prediction, k is 

the predictor level and W., W (k) the respective sample-means 

J s 

of W., W (k) from the 1100 soundings. As seen in Table 1, 

J s 

the arrangement of the predictor levels (denoted by aster- 
isks) takes primary advantage of the specification of the 
W -profile in the lower levels, yet allows for additional 
upper level influences. 

The estimation method presented by WC demonstrated a 

strong multiple correlation between the predicted W. values 

J 

and the indicated saturated mixing ratio as predictands. 
However, in the collection of 1100 sample soundings, strati- 
fication of the regression results of WC by season was not 
considered. Additionally, in application to the computation 
of precipitable water vapor for radiometric (VTPR) calcula- 
tions, island and coastal radiosondes were considered as 
being representative of open ocean areas. 

In this thesis several multiple regression estimation 
formulas for water vapor porfiles are determined based upon 
the Wj and W s (k) data samples from only ten Ocean Station 
Vessels (OSV) distributed over the Atlantic and Pacific. 
However, the OSV data were based only on soundings from the 
month of January in the years 1967 through 1970, and so 
should possibly be considered representative of midwinter 
moisture predictions over the oceans. The regression formu- 
las derived were based upon 1081 dependent data samples and 
then applied for independent test purposes upon a smaller 
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data sample drawn from the same set of OSV-soundings . Eval- 
uation of the various estimation methods derived here was 
conducted in order to arrive at a "best" set of mixing-ratio 
specification equations. This prediction scheme was then 
compared with the WC prediction equations for the purpose of 
revealing the relative specification significance in the 
equations derived here. This evaluation was conducted using 
the smaller independent OSV data sample. Finally, a data- 
sample developed from 160 island and coastal soundings, col- 
lected from 1-28 March 1973 and similar to those used by WC 
(1973), was subjected to the "best" specification system de- 
veloped here and compared with independent testing based on 
the use of the WC specification system. 

A considerable amount of linear interpolation is required 
to adjust sounding data to the 100-level atmosphere, specif- 
ically to the levels 61 through 100. In this thesis and in 

the work done by WC simple linear interpolation on pressure 

"t h ^ 

raised to the 2/7 power as listed in Table 1 is used. It 
is reasonable to assume that the regression coefficients de- 
termined by WC are not independent of one another over the 
40-level system, and so could be reduced to a reasonably 
concise but smaller predictand set, say J=l,...,19. For 
purposes of this thesis, the best-fit multiple regression 
equation in the form of (2) was therefore applied only at 
the 19 levels listed in Table 2. It is to be noted that 
this set of 19 levels is a subset of the 40-level predictand 
set of WC and that it also includes the 11 levels used as 
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Table 2. 19-level atmosphere. The pressure levels P(J) are 
obtained by linear scaling of P 2 / 7 over the interval 
188.45 mb to 1000 mb in 19 increments. Column 3 
indicates the subset of predictor levels. 



Level 


Predictand 
Pressure Level 


W (k) Predictor 
s Levels 


(J) 


P(J) (mb) 


(k) 


1 


188.45 




2 


209.94 


1 


3 


245.37 


2 


4 


271.24 




5 


299.01 


3 


6 


328.75 




7 


360.54 




8 


412.26 


4 


9 


469.08 




10 


509.93 


5 


11 


553.26 




12 


599.16 




13 


647.72 




14 


699.03 


6 


15 


753.17 


7 


16 


810.25 


8 


17 


870.35 


9 


18 


933.56 


10 


19 


1000. 00 


11 
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predictor-levels in their system. It was assumed that re- 
stricting Eq. (2) to 19-predictand levels would result in 
little, if any, loss of resolution of the estimated mixing 
ratio profiles over the lowest 40 WC predictand j-levels. 
Hence, comparison of the regression results deduced here 
will be with those results of WC at the 19-predictand levels 
adopted for use in this study. 

Regression analysis was carried out on the IBM 360 com- 
puter using the BMD03R statistical program (Dixon, 1973). 
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II. DATA PROCESSING 



A. OCEAN STATION VESSEL RADIOSONDES 

The raw data used in this study consisted of a set of 
January radiosondes containing temperature and relative hu- 
midity data in accordance with standard radiosonde reporting 
procedures, i.e. available mandatory level plus any signifi- 
cant levels (Radiosonde Observations Handbook, 1969), from 
the 10 OSVs. Table 3 gives the approximate latitude and 
longitude for each of the OSVs in addition to the sample- 
size contributed by each to the total data set. Also listed 
in Table 3 is a similar breakdown of the sample size contri- 
butions to the total dependent and independent data samples 
(abbreviated as OSV-1, and OSV-2, respectively). The latter 
two data samples will be discussed later in this section. 

Each sounding included in the OSV data samples was re- 
quired to have temperature data up to at least 150 mb and 
relative humidity data to 450 mb. Both mandatory and sig- 
nificant levels were combined to produce temperature and 
relative humidity profiles at pressure levels arranged in 
order from 1000 mb to 150 mb. 

B. REDUCTION OF RADIOSONDE OBSERVATIONS 

Once a radiosonde was arranged in a pressure-decreasing 
order the Goff-Gratch formula in the form (List, 1963) 
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Table 3. Radiosonde-data sample sizes from Ocean Station Vessels, and 
their geographic locations. 



osv 


POSITION 

LATITUDE LONGITUDE 


TOTAL 

OSV 


NUMBER OF 
DEPENDENT 
(OSV-1) 


SOUNDINGS IN SAMPLES 
DATA INDEPENDENT DATA 

(OSV- 2) 


A 


62. 5N 


33. 0W 


194 


97 


97 


B 


56. 5N 


50. 5W 


154 


77 


77 


C 


53. 5N 


35. 5W 


189 


95 


94 


I 


59. ON 


19. 0W 


212 


106 


106 


J 


53. 5N 


19. 5W 


206 


103 


103 


P 


50. ON 


145. 0W 


202 


101 


101 


D 


44. 5N 


41. 0W 


169 


113 


56 


E 


35. ON 


48. 5W 


185 


123 


62 


N 


30. ON 


140 . OW 


205 


136 


69 


V 


38. ON 


165. OW 


194 


130 


64 


TOTALS 




1910 


1081 


829 
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log(e /Po ) = -7 . 90298 ( (T /T )-l )+5 . 02808 log(T /T) 
s s s 

-1.3816 x lO-’ClO 1 1 ' 3 34(1-T / T s)_i) 

+8.1328 x 10 -3 [10 -3 ' 49 1 49( (T s/ T ^ _1) -l] , (3) 

where T = 373.16°K, P 0 = 1013.246 mb, and T = T(P) is the 

o 

temperature at pressure P, was used to determine saturated 

vapor pressures W . From the value of e , W (P) was deter- 

s s s 

mined from 

.622 e 

V P > = P -e S O) 

s 

which in turn yielded the actual mixing ratio W(P) upon appli- 
cation of the relationship between RH (as a fraction) and W 

s 

W(P) = RH x W o (P) (5) 

Following this procedure, simple linear interpolation of 

both T(P) and W(P) relative to the [P(J)] 2 / 7 values of Table 

2 was performed, yielding a temperature profile T(J) up to 

188.45 mb and a mixing ratio profile W(J) to 469.08 mb. 

Equations (3) and (4) were again needed to establish the 

11-predictor values of W (k); in this case T = T(k), P = P(k), 

s 

k=l,...,ll, as shown in Table 2. 

Since mixing ratio (or RH) profiles were usually unavail- 
able in the radiosondes above 450 mb, the well known cubic 
power-law profile of Smith (1966) expressing W in the form 

W = W ..(P/P ,)\ A = 3 (6) 

ref ' ref' ’ 
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was used. Here W ^ and P ^ signify the simultaneous con- 
ditions at the reference level, taken here as P _ = 469.08 mb. 

ref 

Equation (6) was employed to compute "observational" values 
of W at all levels above P re ^ to complete the 19 levels of 
Table 2. 

Completion of the above procedures through (6) yielded, 
for each sounding, an array of 19 values of W and 11 values 
of W . These values were arranged in a final 30-variable 
format [W(1),W(2), . . . ,W(19) ,W g (l) ,W g (2) , . . . ,W g (ll)] . This 
format was found to be most convenient accessing the least 
squares statistical program, BMD03R. 

C. DIVISION OF DATA SAMPLES 

The collection of 1910 OSV soundings was divided into 
two major data samples, OSV-1 consisting of 1081 total 
soundings, and OSV-2, comprised of the remaining 829. OSV- 2 
data were obtained by selecting a random mix of every second 
to every third sounding from the total set of soundings 
available from each ship, with the remainder going into the 
dependent sample OSV-1. Table 3 lists the contributions to 
each sample on a ship by ship basis. OSV-1 was used for de- 
termination of the dependent-test set of regression equa- 
tions, for W( J) specified in terms of the predictor-variables 

W (1),...W (11). OSV-2 was used for independent or verifi- 
s s 

cation tests of regression results developed using OSV-1. 

An additional subdivision of the OSV data was made to 
determine latitudinal effects, if any, on the moisture esti- 
mation regressions. This subdivision was suggested after 



21 



consideration of the distributions of the mean mixing ratios 
W(J) for each ship station in the OSV-1 data sample. 

Table 4 lists the mean W-values by level for the entire 
dependent sample set, the subdivisions of the sample accord- 
ing to 4>_>50N, and 4><50N, as well as by individual ships. It 
seemed plausible to use 50N as the dividing point for the 
latitudinal zones as this arrangement provided reasonably 
large sample sizes to stratify the regression analyses in 
addition to maintaining the apparent stratifications demon- 
strated by the mean distributions. Although Ship Station P 
was located at 50N, its W-profile statistics seemed to be 
described more closely by the means of the northern set as 
its verification statistics suggest (Sec. IV.B.l). 

Verification of the estimation equations derived by the 
latitudinal stratification of the dependent sample set was 
conducted on the stratifications (a) all <f> , (b) $_>50N, (c) 

4><50N as well as for the individual ship-data samples re- 
served for the independent OSV-2 sample. The results of the 
tests on the individual ships are tabulated in Appendix B. 

D. ISLAND AND COASTAL RADIOSONDE SAMPLE 

Data collected under the auspices of the NATIONAL ENVI- 
RONMENTAL SATELLITE SERVICE (NESS) during March, 1973, were 
used for additional testing of the estimation procedures. 

This data set provided a total of 160 typical radiosondes 
containing moisture data from the surface up to at least lev- 
el J=9 of the pressure scale of Table 2 but with temperatures 
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Table 4 (a). Mean W-profiles for OSV-1 data, gm/lOOkg. 

Stratification Individual Ships 

Level All <J> 4> >50 4><50 A B C I J 



> 



2 



w 



Q 



ft 



n- O' ^ sO co <0 cm coin co m* cm -^ co^vO OM- cm 




h- co ^ co O' Mf-o^rx-incMCMoo^^^oococx. 




CM ro sO 00 0 s4- O' CM xO x4- in a> >4 - 0 ^ CM CM ^ 




H^^ninC-CD^lAONst vt sDO 




•h ^cm cm co si* inc- 
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'4‘ooincMcoo'coM*sO>4'»-«r s -xOO'CMCDino\0 




co > 4 * 0 co o' co co 00 ^ in co rx- in O' O ' 00 rx- n4* 
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up to at least level 1. The format of the soundings was 
much the same as that of the OSV data except that the mois- 
ture element was given in terms of the dew point T d rather 
than RH. This difference required the use of Eq. (3) with 
e and T replaced by e and T , , respectively, followed by 
application of Eq. (4) to obtain the W-values. All other 
procedures in this section apply to the NESS data to obtain 

soundings in the form [ W ( 1 ) , W(2 ) , . . . , W( J) , W (1),W (2),... 

s s 

W (11)], where J is at least as high as level 9. Recall that 
s 

W (1),W (2 ),..., W (11) depend simply upon having the temper- 
s s s 

ature profile complete in the form of Table 2 to level 

P(2) = 209.94 mb, and that for purposes of computing W (k), 

s 

the indexing k=l,...,ll denotes the indicated subset of the 
J-values shown in column 3 of Table 2. 



E. VERTICAL WEIGHTING 

With the mixing ratios estimated at most of the 19 levels 
of Table 2, it was possible to compute precipitable water 
vapor for each sounding as a measure of the vertically inte- 
grated mixing ratios. The precipitable water vapor, calcu- 
lated by means of the finite-difference integration, Eq. (7). 



U = 



1<J<9 

Z (W J+1 +W )(P ~P ) 
1000 J=18 j+1 J j+1 _ 

980 2 



(7) 



provided a basis for testing "estimated" precipitable water 
vapor against the corresponding "observed" U-values deduced 
from the W-profiles summed over the same atmospheric levels 
from level J downward. In the case of OSV soundings the 
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summation of Eq. (7) extended up to J=l; however in the case 
of the NESS data, the top level of humidity reports extended 
up to a level J in the range 1<J<9. 
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III. STATISTICAL CONCEPTS AND PROCEDURES 



A. MULTIPLE LINEAR REGRESSION 

This thesis makes use of multiple linear regression 
(MLR) techniques to develop an effective specification 
scheme to estimate mixing ratio values at the 19 previously 
designated pressure levels, using as predictors, the satu- 
rated mixing ratios obtained at the 11 key levels specified 
by WC. In the discussion to follow, the phraseology "de- 
pendent variable" refers to observed mixing ratios to be 
specified in terms of the independent variables when MLR 
techniques are to be applied. The independent variables are 
always the set of saturated mixing ratios computed at the 11 
specified levels in the vertical soundings. 

Assuming that the dependent variable Y depends on k in- 
dependent variables (X^ ,X^ , . . . , X^) , and that there is a sam- 
ple of N observations of the k+1 variables (Y.X^.X^, . . . ,X^) , 
MLR techniques can be applied to find the best-fit multiple 
regression coefficients C.. Y can be estimated by the re- 
gression equation 

k 

Y = A + E C .X . (8) 

J = 1 J 3 

A 

where the superscripted symbol Y denotes "estimation of Y." 

Also it follows, from the requirement that the regres- 
sion plane pass through the sample mean, that 

k 

A = Y - .E C.X.. (9) 

J = 1 3 3 
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The quantities ( ) signify the respective means based on the 
sample size N. 

It should be noted that Eqs. 8 and 9 can be combined to 
produce the form 

k 

Y - Y = .IL C .(X .-X . ) (10) 

J=1 J J J 



where the C.'s are the best-fit regression coefficients de- 

J 

termined by MLR techniques. 

The regression coefficients of (10) are determined by 
solving the system of normal equations given in determinant 
form by 





a 

i i 


a i 2 * ’ * a ik 




a iy 


| C i’C 2 ,.--»C k 


a 2 1 


a 2 2 * - ’ a 2 k 


= 


a 2 y 




a k i 


a kk 




a ky 



where a. . and a. are given by 
ij iy 



a . . 
ij 



N 

= £ (X. -X. )(X . -X . ), 

m=l v lm i /v jm j 



and 



N 



a. = £ (X. -X. ) (Y -Y), 

iy m=l v lm i m 



✓ -* -« \ 



V. J-j. j 



( 12 ) 



(13) 



respectively . 

In Eqs. (12) and (13) it is understood that (Y (m \x[ m \ 

...,x5 m )) constitute the mth sample of Y on X. . For simpli- 
’ k ' 1 

fication the superscript in this sample designation has been 
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written as a subscript in (12) and (13). From Eq. (12), 

a. . = a. . , and therefore, the matrix of {a. .} in (11) is sym- 
^ J J ^ 1 J 

metric. Equation (11) represents a system of k equations in 
the k unknowns which can be solved for a unique set of 
regression coefficients, provided the k-by-k determinant of 
the left hand side of (11) is non-singular. If the multipli- 
cation of the left hand terms of (11) is carried out, one 
can obtain k equations in the form 



a, C 
k i i 



+ a. C 
k2 2 



which can be solved (Crow, 
methods applied to {a^}. 



+ a kk C k “ a ky 
et al., 1955) by matrix 



(14) 

inversion 



B. SIGNIFICANCE TESTS 

1 . On Dependent Regression Equations 

The coefficient of determination of an MLR technique 
is definable as the ratio of the sum of squares due to the 
regression to the total sum of squares. This parameter may 
be shown to be the square of the multiple correlation coeffi- 
cient, R^, 



■r 2 _ Sum Sq. Explained by Reg. 
k M Sum Sq. of Total Deviation 



N 

[ i=i ( V Y)2] 



N _ 

[ i=i ( V Y)2] 



(15) 



and is used as a measure of the efficiency of the regression 
performed by the BMD03R program. The total sum of squares 

N _ N _ 

S.S. total = ^(Y^Y) 2 = i = 1 t(Y i -Y)+(Y i -Y.)] 2 = (N-l)S 2 (16) 
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S.S. total may be expanded 



defines the variance Sy. 



N 



S.S. total = E (Y.-Y) 2 = 
i=l v 1 / 



N N 

. E (Y.-Y) 2 +. E (Y.-Y. ) 2 

1=1 l 1=1 l i' 



( 17 ) 



(a) 



(b) 



which is a combination of the S.S. due to regression (a), 
and the residual S.S. (b). 

If the regression were perfect, term (b) of (17) 
would vanish and R 2 j would be 1; conversely, if Y were totally 
independent of X^, would be 0. Combination of Eqs. (15) 
and (17) shows that the residual sum of squares after appli- 
cation of the regression (8) is 



N 

Res. S.S. = (1 - R 2 ) E (Y.-Y) 2 . 

v M 1=1 l / 



A method of determining the significance of a predic- 
tion equation having k predictors is the F-statistic with k 
and N-k-1 degrees of freedom, defined by 



F (k,N-k-l) 



_ Mean Square Explained by the Prediction 
Mean Square Unexpl. after the Prediction ’ 



(18) 



or 



F(k, N-k-1) = 



V k 



(1-R 2 )/ (N-k-1) 



(19) 



Equation 19 is a direct application of the defini- 
tion (18) of the F-statistic for the appropriate degrees of 
freedom. For the numerator of (19), k is the number of in- 
dependent variables used in specifying the mean square ex- 
plained of (18). In the denominator of (19) the degrees of 
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freedom (df) is equal to N-k-1, corresponding to the total 
sample size, reduced by the 11 known coefficients and the 
value of the sample mean. In the tests on the dependent 
sample k=ll throughout, although N varies from test-sample 
to sample. 

Critical values of the F-statistic, designated F , 
can be found in standard texts on statistics for the appro- 
priate degrees of freedom at specified levels of signifi- 
cance. If the F-statistic of (18) exceeds F. 0 j (k, N-k-1 ) , 
there is less than 1 % probability that the sample came from 
a population in which there is no linear relationship be- 
tween Y and X^. Such a statistical outcome for the regres- 
sion is described as corresponding to a 99% level of 
confidence . 

A third related statistic which is listed as output 
by BMD03R is the standard error of estimate (after applica- 
tion of the regression) is defined by 



(STD. 



err.) 2 = 



S.S. Errors 
(N-k-1) 



Z(Y.-Y. ) 2 

l l ' 

(N-k-1) 



( 20 ) 



Equation (20) shows that the standard error is simply the 

square root of the mean of the errors squared. The standard 

error is related to R, 2 . by 

M 



(STD. ERR. ) 2 



(N-1)S 2 (1-R 2 ) 

(N-k-1) 



( 21 ) 
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The three statistics defined in (15), (18), and (20) 
were considered as comparative measures of the MLR techniques 
of this work and were conveniently provided as a standard 
output from the BMD03R program. Values of R^, STD. ERR. , 
and F are presented in tabular form in Section IV, based 
upon MLR techniques applied to relevant subdivisions of both 
dependent and independent samples. 

2 . On Independent Data for Level by Level Analysis 

Similar test procedures were applied to evaluate es- 
timates of the regression methods from Sec. III.A.l to the 
independent OSV-2 sample. Slight variations in the degrees 
of freedom were required in specifying the efficiency sta- 
tistics defined by Eqs. (15), (l8), and (20) to the inde- 

pendent samples. In the regression analyses applied to the 
dependent sample, 11 predictors and the mean were used in 
arriving at N-12 degrees of freedom in the denominator. 

However once the coefficients were applied to the OSV-2 sam- 
ple, there was essentially only one way to form the estima- 
tions Y, using the prescribed set of coefficients. Therefore, 
in this section the df corresponding to k was taken as 1. 

R^ was also somewhat modified and was termed simply R^p, 
defined as the effective coefficient of determination. The 



R 



EFF 



values were defined by 



R 



EFF 



= 1 - 



S.S. Errors 
S.S. Total 



( 22 ) 
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The F-statistic and the standard error of estimate 
were calculated in a manner similar to Eqs. (19) and (20) 
with k=l, but with N-k-1 replaced by N-l. The resulting ex- 
pression for these statistics 



F 



EFF 



(23) 



and 

Z(Y -Y ) 2 

( STD . ERR . ) 2 = — — ( 24 ) 

were used when conducting significance tests on the inde- 
pendent data set. 

The statistics were collected and tabulated at each 
of the 19 prediction levels. Comparison of their values 
from the various estimation methods provided a means of de- 
termining whether a lat itudinally dependent specification 
scheme over the OSV sounding samples was more effective than 
a composite-latitude specification scheme. After determina- 
tion of a "best" scheme, specification usefulness was com- 
pared in the same manner with the WC prediction scheme. 

(Here WC denotes the Weinreb-Crosby regression system. ) 

3. Independent Data by Precipitable Water Vapor 

Precipitable water vapor values, both actual and 
predicted, calculated by Eq. (7) were subjected to a simple 
linear regression of U on U by BMD03R, in order to obtain a 
vertically weighted comparison of the various specifications. 
The three statistics described in Section III.B.l were used 
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to judge the effectiveness of the estimates by the various 
methods . 
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IV. STATISTICAL RESULTS 



A . DEPENDENT SAMPLE 

This phase of the research consisted of collecting, 
tabulating, and evaluating the results of three MLR tests 
carried out by BMD03R on the dependent OSV-1 sample of sound- 
ings. The first MLR provided the 12 coefficients and the 
statistics for each of the 19 prediction equations to be 
applied over the composite sample of ship-soundings. The 
second and third MLR's provided similar sets of coefficients 
and statistics for the latitudinal divisions of the soundings 
for <(>_>50N and <})<50N, respectively. These three sets of re- 
gression equations, employed in the form of Eq. (8), will be 
referred to as: 

1) COEF-1 for the one over-all set, regardless of <j> , 

2) C0EF-2N,S using 

a) The coefficient set C0EF-2N when <}>>50N, and 

b) C0EF-2S when <f><50N, 

3) The composite set of COEF-2 without the N or S 
designation . 

The latter set implies a composite estimation by one or 
other of the two components of the applicable equations in 
the appropriate latitudinal zones. The sample size is the 
same as for COEF-1. 

Table 5 lists the values of R^ , the standard error of 
estimate, and F as obtained from BMD03R for each MLR experi- 
ment at each of the 19 levels listed in Table 2. The most 
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Table 5. Summary of statistics obtained during development of MLR coefficient sets COEF-1, C0EF-2N, 
and C0EF-2S by BMD03R from dependent OSV-1 sounding samples. 
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stringent critical F-value, F (k,N-k-l), for the entire dis- 

o 

play at the 1% significance level occurs for the degrees of 
freedom indicated by 

F (11,490) = 2.98, 
c 

which occurs in the southern subdivision of OSV-1. As can 
be seen from Table 5, all F-values exceed this critical val- 
ue by a rather impressive amount. The inference was drawn 
that each estimation procedure so determined was based on a 

strong linear dependence of W on W (k). Additionally, com- 

s 

parison of the standard errors of Table 5 with the corre- 
sponding OSV-1 sample standard deviations (listed in Table 
4) indicates a significant reduction of the variances by the 
regression estimates. 

An example of a regression equation at level J=19 (cor- 
responding to 1000 mb) is given by 

W(19) = 2.1082+0.4979 W g (l)-0.5900 W (2)-1.1938 W (3) 

-0.1600 W (4)+0. 2194 W (5)-0.2095 W (6)-0.1151 W (7) 
s s s s 

+0.2814 W (8)+0. 0815 W (9)+0.9076 W (10)-0.0313 W (11) (25) 

s s s s 

Equation (25) is one prediction equation of the set of 19 
equations developed from COEF-1. A listing of the three co- 
efficient sets comprising COEF-1, C0EF-2(N,S) as well as the 
WC coefficient set is contained in Appendix A. It should be 
noted that the first term on the right side of (25) results 
from Eq. (9) with all the sample means known. 
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Examination of Table 5 does not reveal convincing evi- 
dence with which to determine a "best" prediction method be- 
tween COEF-1 and COEF-2 as all F-values demonstrate at least 
a 1% significance level for dependent W-specif icat ions . For 
purposes of this study, a highly significant improvement in 
specifications by COEF-2 over COEF-1 would be required. This 
requirement was deemed essential to override the added com- 
puter storage and procedures required to implement a double 
set of coefficients. Clearly, Table 5 does not offer this 
evidence . 

B. INDEPENDENT SHIP SAMPLE 

1 . An Illustration of the Usefulness of Lat itudinally- 

Stratified Coefficient-Set 

Recall that in Section II, the question of a latitu- 
dinal stratification of the water vapor prediction equations 
was raised. The reason for considering this approach is 
well illustrated by the level-by-level results for Ship P 
presented in Table 6 (based here upon independent test data). 
The results of Table 6, e.g. larger values of R EFF resulting 
from the C0EF-2N regression scheme, indicated stronger sig- 
nificance level-by-level than those of either COEF-1 or 
C0EF-2S coefficient sets when applied to the same data. 

These results, as much as any individual ship results, indi- 
cated the need to consider, initially, a lat itudinally- 
stratified set of predictors, with the boundary set at b^SON. 
In addition, COEF-1 regressions consider the alternate hy- 
pothesis, namely, the case of no latitudinal stratification. 
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Table 6. Verification statistics resulting from application of COEF-1, C0EF-2N, and 
C0EF-2S coefficients to the 101 OSV-2 soundings from Ship P. 
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V 




Note that like results for the remaining nine ships 
of the OSV-2 sample are listed in Appendix B in the order 
specified by Table 3. 

2 . For the Determination of a "Best" Set of Coefficients 
a. Level-by-Level 

Y/ater vapor estimates were made by COEF-1 and 
C0EF-2N.S for the entire independent OSV-2 sample. The re- 
sults are to be tested level-by-level as follows: 

i. COEF-1 vs. C0EF-2N for latitudes <f>j>50N 

ii. COEF-1 vs. C0EF-2S for latitudes <(><50N 

iii . COEF-1 vs. COEF-2 for all $ 

The statistical results for the three comparisons are listed 
in Tables 6 through 8, respectively. As in Section IV. A the 

F value for a 1% significance level was taken for the most 
c 

stringent case, i.e. case (ii) was 
F (1,250) = 6.75 . 

The summary of F-statistics in Table 7 for <f>>50N 
shows strong predictive capability at all levels in both 
COEF-1 and C0EF-2N, with the latter demonstrating slightly 
higher verification (F-value) statistics. The results of 
C0EF-2S were not as strong when compared to COEF-1. Table 8 
show's slightly better verifications with COEF-1 from level 1 
through 11, then a reversal for the remaining eight levels. 
Note that at level 9, for the C0EF-2S tabulation, the value 

R EFF =0, but als ° F<0 * The R EFF -Value ’ aS defined b y Eq ‘ (22 ^ 
takes on negative values whenever the mean square error of 
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Table 7. Statistical comparison of independent data verification by COEF-1 and C0EF-2N for <}>>50. 
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Table 8. Statistical comparison of independent data verification by COEF-1 and C0EF-2S for (|)<50. 




41 



the estimation scheme exceeds the sample variance. For this 
case, the sample-mean must be taken as the best prediction 



scheme, and is reset to 0 in obtaining a realistic 

value of Rppp- However the F-value was left as calculated 
by Eq. (23), with k=l and N-l degrees of freedom. The sta- 
tistical comparison of COEF-1 and COEF-2, listed in Table 9, 
for the entire OSV-2 independent sample was a composite of 
the results of Tables 7 and 8. Both estimation methods 
demonstrated strong prediction qualities, as evidenced by 
the large F-values. As in the regression analyses in Section 
IV. A, the comparison of the performance of COEF-1 and COEF-2, 
Table 9, is quite close over most of the 19 levels. Hence, 
there was no conclusive evidence supporting a decision to 
select a ''best" set of coefficients from either COEF-1 or 
COEF-2. This decision was deferred to the outcome of the 
regression results relating observed to estimated values of 
precipitable water vapor. 

b. Precipitable Water 

Precipitable water calculations (U) were made on 
the OSV-2 sample soundings according to Eq. (7). Similar 
application of (7) to the 19-level OSV-2 mixing ratio speci- 
fications generated by COEF-1 and the composite set of 



COEF-2 equations produced two modes of estimated U-values 
for each of the 829 soundings, depending upon the latitude 
of the sounding. 

The simple linear regression analyses were per- 
formed on each set of U, U using BMD03R. Each regression, 
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Table 9. Statistical comparison of independent data verification by COEF-1 and composite COEF-2 
for all cj). 
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identical to the MLR with only one independent variable, was 
performed with U as the dependent variable and the appro- 

A 

priate U as the independent variable. These regressions 
provided significance measurements of the 19-level predic- 
tions by a single vertically-scaled F-statistic with 1 and 
N-2 degrees of freedom. 

The statistical results from the precipitable 
water vapor analyses are listed in Table 10 and provided a 
final comparison between the COEF-1 and C0EF-2 specification 
systems. The critical F-value, 

F = 7.88 
c 

corresponding to the sample of N=251 soundings (<f><50N) pro- 
vided the most stringent F-statistic significance require- 
ments namely at the 1% level. As seen in Table 10, all 
calculated F-values greatly exceeded the critical value, 
thus indicating, to at least the 1% level of significance, 
that none of the specifications of precipitable water could 
have occurred by chance. 

The precipitable water results were analogous to 
those in Tables 7 through 9, and as seen in Table 10, the 
single COEF-1 set of equations yielded the most effective 
specifications. Considering the conditions set down in Sec- 
tion IV.B.l requiring considerably better specification by 
COEF-2 over COEF-1, and with the results of Sections IV. A 
and IV. B, the decision was made at this point to eliminate 
from further consideration the model of a lat itudinally 
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differentiated set of coefficients. Thus the remainder of 



this research incorporated only the composite COEF-1 set of 
equations for comparison with the WC specification method. 

A 

Table 10. Precipitable water vapor regression statistics of U on U by 
COEF-1 and by the appropriate C0EF-2 scheme applied to inde- 
pendent 0SV-2 sample. 







COEF-1 


C0EF-2 (N, S , Composite) 


AREA 


MEAN STD 

DEV 


STD R F 

ERR R F 


STD R F 

ERR R F 


<J>>50 

(N=578) 


.998 .408 


.2344 .8191 1170.1 


.2418 .8054 1063.8 


4><50 

(N=251) 


1.73 .674 


.3847 .8218 517.9 


.3902 .8161 496.4 


All 4> 
(N=829) 


1.22 .606 


.2881 .8800 2838.6 


.2947 .8741 2677.2 



Operational implementation of a single set of 
coefficients over a more complex set eliminated the require- 
ment of additional computer core allocation for storage of a 
second set of coefficients and the necessary logic state- 
ments to enact the stratification differentiation. Clearly, 
the current results do not warrant use of the more complex 
COEF-2 prediction scheme. 

3 . Comparison With the WC Predictors 
a. Level-by-Level 

In this section comparison of estimates produced 
by the "best" COEF-1 coefficient set was made with those es- 
timates from the WC predictions. Comparison was conducted 
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in accordance with the precedures discussed in Section 

IV.B.l.a. As in previous comparisons, the most stringent 

conditions were adhered to for determination of an F -value, 

c 

in this case 

F (1,250) = 6.75 . 

Statistical comparison of the two specification methods are 
listed in Tables 11, 12, and 13 for the lat itudinally strati 
fied subsamples (0SV-2N, 0SV-2S) and for the total OSV-2 
sample, respectively. 

Without exception, the COEF-1 specifications 
were statistically superior to those produced by the WC co- 
efficients at all 19 J-levels in each of the three tables. 
The WC specifications do demonstrate considerable predicta- 
bility in the lower levels of the atmosphere, however, the 
F-values indicate poor results in the highest two or three 
levels. Some possible explanations will be offered in Sec- 
tion V for the breakdown of the upper level specifications 
by the WC equations. It is clear from the results listed in 
Tables 11 through 13 that the COEF-1 specifications are su- 
perior to the WC predictions for the level-by-level analysis 
b. Precipitable Water 

A final comparison of the WC and COEF-1 specifi- 
cation methods was made using precipitable water vapor cal- 
culations as in Section IV.B.l.b. Table 14 lists the 
statistics resulting from the vertically-scaled evaluation 
of the two water vapor estimates for each of the OSV-2 
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Table 11. Statistical comparison of independent data verification by COEF-1 and WC for <f>>50. 
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Table 12. Statistical comparison of independent data verification by COEF-1 and WC for 4><50. 
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Table 13. Statistical comparison of independent data verification by COEF-1 and WC for all <J> 
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453 



stratifications. The statistics were generated from two 
simple linear regressions by BMD03R of U on each of the U- 
estimates. The last line of Table 14 indicates a definite 
' superiority of the COEF-1 equations over the WC equations 
for the independent OSV-2 data sample. 

Table 14. Precipitable water vapor results by COEF-1 and modified WC 
coefficients applied to independent OSV-2 sample. 







COEF-1 


WC 


AREA 


MEAN 

* DEV 


STD R F 

ERR R F 


STD R F 

ERR R F 


<}>>50 

(N=578) 


.998 .408 


.2344 .8191 1170.1 


.2658 .7587 781.1 


<(><50 

(N=251) 


1.73 .674 


.3847 .8218 517.9 


.5324 .7297 283.5 


All <p 
(N=829) 


1.22 .606 


.2881 .8800 2838.6 


.3393 .8289 1816.1 



C. ISLAND AND COASTAL (NESS) DATA 
1 . Level -by- Level 

The 160 NESS island and coastal March soundings pro- 
vided a test sample that focussed, to a degree, on evaluation 
of the possible significance of winter effects on W-profile. 
Table 15 presents the statistics obtained by the first ef- 
forts at estimating the W-profiles by COEF-1 and WC coeffi- 
cients. In making these first W-specif icat ions all 160 
soundings (with observed W-values up to at least 510 mb) 
were accepted as valid data. The results presented in Table 
15 were at best discouraging for both sets of coefficients. 
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Table 15. Statistical comparison of COEF-1 and WC coefficients applied to the 160 NESS soundings. 
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The negative F-values, starting at level 13 for the COEF-1 
specifications and at level 9 for those by the WC method, 
indicated considerably greater effectiveness would be achieved 
by making predictions with the sample mean values of the mix- 
ing ratios. 

It was noted that a number of the NESS soundings 
should not have been considered reliable after reviewing the 
actual W-profiles. Unreliable soundings were identified 
when the reported moisture profiles indicated RH<20% in the 
lower levels of the soundings and persisted with increasing 
ascent. It was felt that the low moisture reports were due 
mostly to a failure of the moisture element on the radio- 
sonde instrument. 

If a sounding had RH considerably higher than 20% at 
all levels, its moisture profile was considered reliable and 
the sounding was retained. In some cases the relative hu- 
midity decreased to values lower than 20% in the mid-tropo- 
sphere but then recovered to exceed 20% at higher levels. 

Such soundings were also retained in a revised test set. If 
the RH persisted at less than 20% throughout most of the 
troposphere, including the upper regions where motor-boating 
is not infrequent, the sounding was considered unreliable 
and discarded (see Table 16). 
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Table 16. Examples of three types of soundings categorized 
by RII profiles extending up to level 10 and con- 
sidered for possible inclusion in the NESS sample. 



Level 


Reliable 




(% RH) 
Recovered 


Unreliable 


10 


72 




71 


15 


11 


67 




63 


16 


12 


62 




23 


16 


13 


45 




11 


17 


14 


42 




9 


18 


15 


43 




11 


21 


16 


42 




16 


20 


17 


37 




31 


18 


18 


57 




29 


67 


19 


79 




38 


76 


After 


screening 


the 


NESS soundings 


using the method 


of Table 16, 


a total of 


124 


soundings remained for testing. 



Table 17 lists the statistical results of the level-by-level 
W-specif ications by the COEF-1 and WC equation sets. Since 
the number of soundings with observed W-values drop off 
sharply above level 10, the critical F-value was arbitrarily 
taken for N=70 occurring at level 7 and 

F =7.01. 
c 

This value of F is considerably more stringent than neces- 

c 

sary for specifications at levels below level J=7 , but as 
Table 17 demonstrates, this is a reasonable restriction. 



53 



Table 17. Statistical comparison of COEF-1 and WC coefficients applied to the 124 NESS soundings. 
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Examination of the F-values of Table 17 indicated 
that the WC coefficients performed more effectively in the 
lower levels of the atmosphere. Above level 16, the COEF-1 
set of equations produced statistically superior estimates 
up to level 4, where both estimation procedures produce neg- 
ative and F-values. The small number of soundings con- 

taining data in the last four levels was considered responsible 
for these unusual statistical outcomes and it was impossible 
to draw any useful conclusion regarding the comparative 
upper-level specificity by the two estimation systems. 

2 . Precipitable Water Vapor 

Precipitable water vapor integrations were performed 
using Eq. (7) on each of the 124 NESS soundings, incorpo- 
rating all available J-levels for which mixing ratio values 

A 

existed, J<9. Additionally, U-values were calculated using 

A 

only the corresponding levels for the associated U-integra- 
tions up to level J. 

The statistical results, as produced by a BMD03R 

A 

simple linear regression of U on each of the two U's, are 
displayed in Table 18. 

Table 18. Statistical results from BMD03R simple linear 

regressions on precipitable water vapor calcula- 
tions obtained from the 124 NESS observed and 
COEF-1, WC-estimated mixing rat io- values . 









COEF-1 






WC 




MEAN 


STD 

DEV 


STD 

ERR 


R 


F 


STD 

ERR 


R 


F 


1.44 


1.1 


.4292 


. 9218 


690.1 


.4850 


. 8990 


514.0 
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The critical F-value 



F ( 1,122) = 6.84 

was compared to the calculated F-values of Table 18 to eval- 
uate the significance of the vertically scaled mixing ratio 
specifications. It was readily apparent that the U-estimates 
both by the COEF-1 and WC systems demonstrated high statis- 
tical significance and that the COEF-1 estimates indicated 
somewhat higher predictability. 

Reviewing the results of Table 18 along with those of 
Table 17, it became obvious that the near-surface mixing 
ratio values dominate the precipitable water vapor integra- 
tion, and that the upper level breakdown of the WC specifi- 
cations, as evidenced by Table 17, was masked by this 
dominance. Note that the W-values of the dependent and in- 
dependent OSV data sets were determined by Smith's power law, 
Eq. (6), for levels above J=9. On the other hand, the W- 
values for the NESS data set were taken solely as indicated 
by the radiosonde instrument, as were the data used by WC in 
determining their regression coefficients. Application of 
the power law to the OSV data tended to smooth the W-profile 
at higher levels and yielded a degree of consistency in the 
OSV W-profiles in the highest 5 prediction levels (P_<300 mb), 
giving small W-values over the OSV's at levels higher than 
J=5 as compared with those resulting from the NESS sounding 
samples (Table 17). It is felt that the high values of W in 
the NESS data could be attributed largely to inaccuracies in 
the radiosonde observational data at the higher levels. 
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V. CONCLUSIONS 



This thesis has reinforced the findings of Weinreb and 
Crosby (WC) in providing supporting evidence that mixing 
ratios can be described by a linear combination of a desig- 
nated 11-element temperature T(P)-prof ile , which is equiva- 
lent to a profile of W g (P). Further, it has been shown that 
any dependence of the relationship between temperature and 
mixing ratio on latitude is insignificant, and that one 
single set of specification equations can be employed over 
all latitudes. 

It was evident from the statistical results presented in 
Section IV. B. 2 that the performance of the coefficient set 
(COEF-1) was superior to that of the WC coefficient set over 
the Ocean Station Vessels. However, the causes of the im- 
proved specifications were not regarded as a result of a 
land-mass (coastal) or seasonal effect. Rather, inaccuracies 
inherent in radiosonde moisture data in the higher level 
atmospheric structure may have affected the definition of 
the WC coefficients, which were derived by moisture data of 
accuracy comparable to that of the NESS data. This reason- 
ing is supported by the independent tests of the COEF-1 and 
WC specifications on the NESS data. 

Efforts to produce an accurate prediction scheme by MLR 
techniques, as demonstrated by this work, appear extremely 
promising. Effective prediction procedures seem not to be 
associated with seasonal, latitudinal or land-mass (coastal) 
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effects. Rather, effective prediction appears to depend 
upon an accurate data base of sufficiently good quality from 
which to define accurate coefficient sets over the depth 
1000-200 mb. The results of this research indicate that the 
OSV data provide such a data base, and that mixing ratio es- 
timates by a set of equations of the form derived here pro- 
vide a representative profile of W(P) up to approximately 
200 mb. The high level of predictability of the precipitable 
water vapor resulting from the vertical integration of mixing 
ratio estimates provided substantial evidence supporting the 
conclusion that the OSV data define highly effective regres- 
sion equations. 
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APPENDIX A 



This Appendix contains a listing of the three sets of 
regression coefficients COEF-1, C0EF-2N, and C0EF-2S deter- 
mined by BMD03R from the OSV-1 ship sample. A fourth list- 
ing contains the selection of coefficients extracted from 
the WC set described in Section I. 

Each of the four sets of coefficients was applied accord- 
ing to the following matrix operation to produce the appro- 
priate mixing ratio estimates at level J: 



1 



| A J’ C J1 ,C J2 



9 * 9 ^ 



Jll 




J 



WsUD 



J = 1, . . . , 19 
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APPENDIX A-l: The coefficient-set (COEF-1) derived from the dependent set of Janaury ship-soundings, over 

all $ , and covering estimations at 19 mixing ratio prediction levels. 
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APPENDIX A-3: The coefficient set (C0EF-2S) derived from the dependent set of January ship-soundings, for 

<(><50N, and covering estimations at 19 mixing ratio prediction levels. 
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APPENDIX A-4: The set of coefficients extracted for the WC coefficients and used for comparative testing. 
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APPENDIX B 



Appendix B lists the level-by-level statistical results 
obtained by application of coefficient sets COEF-1, C0EF-2N, 
and C0EF-2S to the individual OSV-2 soundings. The listings 
are arranged in the order specified by Table 3, only that 
Ship Station P is not included here (see Table 6). 
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APPENDIX B-l: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship A (N=97). 
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APPENDIX B-2: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship B (N=77). 
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APPENDIX B-3: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

and C0EF-2S to OSV-2 soundings from Ship C (N=94) . 
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APPENDIX B-4: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship I (N=106) . 
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APPENDIX B-5: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship J (N=103) . 
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APPENDIX B-6: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship D (N=56) . 
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APPENDIX B-7 : Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship E (N=62) . 
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APPENDIX B-8: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship N (N=69). 
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APPENDIX B-9: Verification statistics resulting from application of COEF-1, C0EF-2N, and 

C0EF-2S to OSV-2 soundings from Ship V (N=64). 
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